Pairing Interaction and Self-Consistent Densities in Neutron-Rich Nuclei 



j3,6 



o 
o 



o 
o 
m 
o 

o 



X 



J. Dobaczewski, '"^ W. Nazarewicz, ' and P.-G. Rcinhard 
^Institute of Theoretical Physics, Warsaw University, ul. Hoza 69, PL-00681, Warsaw, Poland 
Joint Institute for Heavy Ion Research, Oak Ridge National Laboratory, P. O. Box 2008, 

Oak Ridge, Tennessee 37831 
^Institute for Nuclear Theory, University of Washington, Seattle, Washington 98195 
'^Department of Physics and Astronomy, University of Tennessee, Knoxville, Tennessee 37996 
^Physics Division, Oak Ridge National Laboratory, P.O. Box 2008, Oak Ridge, Tennessee 37831 
^Institut fiir Theoretische Physik, Universitdt Erlangen, Staudtstr. 7, D-91058 Erlangen, Germany 

Particle and pairing densities in spherical even-even neutron-rich nuclei are studied within the 
Skyrme-Hartree-Fock-Bogoliubov approach with the density-dependent pairing interaction. The 
influence of the density dependence of the pairing interaction on asymptotic properties of nucleonic 
distributions are analyzed. It is demonstrated that the size of the neutron halo dramatically depends 
on the behavior of the pairing interaction at low density. 

PACS number(s): 21.10.Dr, 21.10.Ft, 21.10.Gv, 21.60. Jz 



I. INTRODUCTION 

One of the main avenues addressed by radioactive ion 
beams is the evolution of nuclear structure as a func- 
tion of neutron-to-proton asymmetry. From a theoreti- 
cal point of view, exotic nuclei far from stability offer a 
unique test of those components of effective interactions 
that depend on the isospin degrees of freedom. Since ef- 
fective interactions in heavy nuclei have been adjusted 
to stable nuclei and to selected properties of infinite nu- 
clear matter, it is by no means obvious that the isotopic 
trends far from stability, predicted by commonly used 
effective interactions, are correct. Such investigations 
reach inevitably into the regime of weakly bound nuclei, 
and these are much more difhcult to treat theoretically 
than well-bound systems Their Fermi energy lies 
very close to zero, and the decay channels must be taken 
into account explicitly. Correlations due to pairing, core 
polarization, and clustering become crucial. In a drip- 
line system, the pairing interaction and the presence of 
skin excitations (soft modes) could require going beyond 
the picture of a nucleon moving in a single-particle orbit 
|^|-^ . The positive aspect of this situation is that drip- 
line nuclei are also critical probes to understand and to 
further develop the nuclear pairing force. 

Surprisingly, rather little is known about the basic 
properties of the pairing force. Up to now, the micro- 
scopic theory of the pairing interaction has only seldom 
been applied in realistic calculations for finite nuclei. 
A "first-principle" derivation of the pairing interaction 
from the bare NN force still encounters many problems 
such as, e.g., treatment of core polarization Hence, 
phenomenological pairing interactions are usually intro- 
duced. In most older nuclear structure calculations, the 
pairing Hamiltonian has been approximated by the state- 
independent seniority pairing force, or schematic multi- 
pole pairing interaction Such oversimplified forces, 
usually treated by means of the BCS approximation, per- 
form remarkably well when applied to nuclei in the neigh- 



borhood of the stability valley, but they are inappropriate 
(and formally wrong) when extrapolating far from stabil- 
ity. The self-consistent mean-field models have mean- 
while reached such a high level of precision that one 
needed to improve on the pairing part of the model. 
The presently most up-to-date models employ local par- 
ing forces parametrized as contact interactions |p^ , pd| ,p| . 
More flexible forms attach a density dependence to the 
pairing strength jl^-Q . There exist even more elaborate 
forms which include also gradient terms [ p"5p^ . It is clear 
that some density dependence is needed. Nuclear mat- 
ter calculations and experimental data on isotope shifts 
strongly suggest that pairing is a surface phenomenon 
and that the pairing interaction should be maximal in 
the surface region. It is not so obvious, however, how 
the actual density dependence should be parametrized 
in detail [Q. This is why neutron-rich nuclei play such 
an important role in this discussion. Indeed, because of 
strong surface effects, the properties of these nuclei are 
sensitive to the density dependence of pairing. Last, but 
not least, a good understanding of pairing is important 
for astrophysical applications where exotic nuclei play a 
crucial role and where the nuclear pairing interaction is 
also important for theories of superfluidity in neutron 
stars [0. 

The main objective of this work is to investigate the 
role of density dependence of pairing interaction on prop- 
erties of neutron-rich nuclei. The material contained in 
this study is organized as follows. The definitions per- 
taining to self-consistent densities and mean fields and a 
specification of interactions are given in Sec. 0. The anal- 
ysis of HFB densities and mean-field potentials obtained 
in different pairing models is given in Sec.|l|. The results 
of self-consistent calculations for pairing gaps, separa- 
tion energies, and halos are discussed in Sec. IV. Finally, 
Sec. contains the main conclusions of this work. 
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II. HARTREE-FOCK-BOGOLIUBOV DENSITIES 
AND MEAN FIELDS 

This section contains a very brief description of the 
HFB formahsm in the coordinate representation. Since 
the method is standard, our discussion is hmited to the 
most essential definitions and references. For more de- 
tails, we would like to refer the reader to the previous 
work 

The HFB approach is a variational method which uses 
nonrelativistic independent-quasiparticle states as trial 
wave functions ^ . The total binding energy of a nucleus 
is obtained self-consistently from the energy functional 



£ — ^kin + ^Sk + £c + £\ 



pair; 



(1) 



where fkin is the kinetic energy functional, Esk is the 
Skyrme functional, £c is the Coulomb energy (including 
the exchange term in the Slater approximation), and fpair 
is the pairing energy. 

By minimizingthe functional (|l|), one arrives at the 
HFB equation |l^Jl^,|^ for the two-component single- 
quasiparticle HFB wave function {01,02}: 



J ^\ h{ra,r'a') -h{ra,r'a') ) \ ME,r'a') 



E + X 
E- X 



h{E,ra) 
h{E,ra) 



(2) 



where A is the Fermi energy and h and h are the particle- 
hole (p-h) and particle-particle (p-p) mean-field Hamil- 
tonians. 

Properties of the HFB equation in the spatial coor- 
dinates, Eq. (H), and the asymptotic properties of HFB 
wave functions and density distributions have been ana- 
lyzed in Refs. |p^ , p^ ,p| . In particular, it has been shown 
that the spectrum of quasi-particle energies E is contin- 
uous for |£'|>— A and discrete for |i?|<— A. Since for A<0 
(bound system) the lower components (t)2{E,ra) are lo- 
calized functions of r, the particle and pairing density 
matrices, 

p{ra,r'a')= ^ MEn,ra)r2{En,r'a') 

0<£;„<-A 

/oc 
dn{E)(b2{E,ra)4>*2{E,r'<j'), (3a) 

p(rc7,rV)-- J2 4>2{En,ra)^l{En,r'(T') 

0<£;„<-A 



dn{E)4>2{E,r(j)c^l{E,r'(T'), (3b) 



are localized as well. For the case of a discretized contin- 
uum, the integral over the energy reduces to a discrete 
sum 



The p-h mean-field Hamiltonian can be expressed 
through the kinetic energy and the p-h mean-field po- 
tential: 



h{ra, r'a') = T{r, r')6„^,, + T{ra, r'a'), (4) 



where 



r{ra,r'a') = J dVad^r'a ^ V^(ra, ratJa; rV, r'acr^) 



PirW2,r2(T2) 



(5) 



and V is the two-body density-dependent effective p-h 
interaction (in our case it is a contact Skyrme interac- 
tion). The pairing mean-field potential can be expressed 
through the p-p density 

h{r<7,r'a') = j A^r'^A^r'^ ^ 2ct'(t^ x 

Vpair(rcr,r', -ct'; r'^CT'i , r'2 , -o-2)/5(ricri , r^cr2) > (6) 

where Vpair is the two-body density-dependent effective 
p-p interaction. (If one adopts the philosophy of the en- 
ergy density functional theory, V and V^air can be dif- 
ferent since they result from different variations of the 
energy functional.) 

The local HFB densities discussed in this work. 



Pir) = ^p(rCT,rcr), 

<7 

Pir) = ^p{ra,ra), 



(7a) 
(7b) 



have a simple physical interpretation ||^. Namely, p{r) 
represents the probability density of finding a particle at 
a given point. On the other hand, |/5p(r) is the proba- 
bility of finding a pair of nucleons in excess of the prob- 
ability of finding two uncorrelated nucleons. 

In this work, in the p-h channel the SLy4 Skyrme 
parametrization pl| ] has been used. This force performs 
well for the total energies, radii, and moments, and it is 
also reliable when it comes to predictions of long isotopic 
sequences pit . 

In the particle-particle (p-p) channel, we employ the 
density-dependent contact interaction. As discussed in a 
number of papers, see e.g. Refs. 0,0,1 , the presence of 
the density dependence in the pairing channel has conse- 
quences for the spatial properties of pairing densities and 
fields. The commonly used density-independent contact 
delta interaction, Vp^ij.{r,r') — V()S{r — r'), leads to vol- 
ume pairing. A simple modification of that force is the 
density-dependent delta interaction (DDDI) |l^-fi4| 



<''h.(^,0 = /pair(r)<5(r-r'), 
where the pairing-strength factor is 

/pah(r)=Fo{l-[pis(r)/pe]"} 



(8) 



(9) 
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and Vq, Pc, and a are constants. In Eq. (||) 
Pis(''') stands for the isoscalar single-particle density 
Pis{f)—Pn{''')+Pp{r). If Pc is chosen such that it is close 
to the saturation density, Pc^pisif = 0), both the re- 
sulting pair density and the pairing potential h{r) are 
small in the nuclear interior, and the pairing field be- 
comes surface-peaked. By varying the magnitude of the 
density-dependent term, the transition from volume pair- 
ing to surface pairing can be probed. A similar form 
of DDDI, also containing the density gradient term, has 
been used in Refs. 15 lq|. 

Apart from rendering the pairing weak in the interior, 
the specific functional dependence on pis used in Eq. (j^) 
is not motivated by any compelling theoretical arguments 
or calculations. In particular, values of power a were 
chosen ad hoc to be either equal to 1 (based on simplic- 
ity), see e.g. Refs. |2^j2^, or equal to the power 7 of 
the Skyrme-force density dependence in the p-h channel 
[pl],p[. The dependence of results on a was, in fact, never 
studied. In the present paper, we perform such an anal- 
ysis by choosing four values of a=l, 1/2, 1/3, and 1/6 
that cover the range of values of 7 used typically for the 
Skyrme forces. 

Calculations which are based on the contact force, such 
as that of Eq. (||) , require a finite space of states in the 
p-p channel. Usually, one takes a limited configuration 
space determined by a cut-off in the single-particle en- 
ergy or in the single-quasiparticle energy. In this work, 
we made a cut-off with respect to "equivalent-spectrum" 
single-particle energies obtained from HFB quasiparti- 
cle energies and occupation coefficients, as defined in 
Ref. [^. All the quasiparticle states with j < jmax=21/2 
and equivalent energies up to emax^GOMeV were consid- 
ered, and the HFB equations were solved in the spherical 
box of i?box=30fm. 

This procedure differs from the prescription applied in 
Ref. [|l^ where a (j£)-dependent cut-off was used. There, 
the quasi-particle states with quasiparticle energies up to 
the depth of the effective potential Dji=—um\Ucii{r\it) 

r 

were considered, and in addition, at least one quasi- 
particle state was considered in each block. The 
combination of these two rules, together with the value 
of i?box=20fm used there, ensured that low-lying high- 
j resonances were always included in the phase space. 
However, systematic calculations across the complete ta- 
ble of nuclides revealed several pathological cases where 
a higher-lying resonance was missing from the included 
phase space. This could be especially true for some 
near-the-barrier j—l—1/2 resonances, for which the ef- 
fective potentials Ucs may have small depths. Since in 
the present study we aim at a better description of the 
continuum phase space, we increased the size of the box 
to i?box=30fm and changed the cut-off prescription so as 
to include sufficiently many states and to never miss a 
resonance. (For more discussion pertaining to the energy 
cut-off problem, see Appendix B of Ref. Q.) 

For Pr we took the standard value of 0.16 fm^"^, and 



the strength Vq of DDDI was adjusted according to the 
prescription given in Ref. |ll]| , i.e., so as to obtain in 
each case the value of 1.256 McV for the average neu- 
tron gap in i20sn^ poj. q,= i^ 1/2, 1/3, and 1/6 the 
adjusted values are yo=-520.5, -787.7, -1041.3, and 
— 1772.5Mevfm~^, respectively. The resulting pairing- 
strength factors (^) are shown in Fig. ^ as functions of 
density p for the four values of the exponent a. It is seen 
that for yO>0.04fm~'^ the pairing-strength factor /pair is 
almost independent of the power a. At low densities, 
however, the pairing interaction becomes strongly depen- 
dent on a and very attractive at p— >0. The pattern shown 
in Fig. |l| indicates that pairing forces characterized by 
small values of a should give rise to pair fields peaked at, 
or even beyond, the nuclear surface (halo region) where 
the nucleonic density is low. 



III. LOCAL HFB DENSITIES AND MEAN-FIELD 
POTENTIALS IN NEUTRON-RICH NUCLEI 

Local p-h and p-p densities are basic elements of 
the HFB-Skyrme theory: they determine self-consistent 
fields; hence the static properties of the nucleus such 
as the binding energy, radius, and shape. The particle 
and pairing local HFB-|-SLy4 neutron densities Pn(r) and 
Pri (?■) calculated for several values of a are displayed in 
Fig. I for isosn. 

With decreasing a, the p-p density p„(r) develops a 
long tail extending towards large distances. This is a di- 
rect consequence of the attractiveness of DDDI at low 
densities when a is small. While in the nuclear inte- 
rior, the p-h density pn{r) depends extremely weakly on 
the actual form of the pairing interaction; its asymptotic 
values are significantly increased when a gets small (see 
inset). Moreover, one observes a clear development of a 
halo structure, i.e., a smooth exponential decrease, that 
for a=l starts at r~6fm, for small a is interrupted at 
r~9 fm, and replaced by a significantly slower decrease 
of the density. 

This is a direct consequence of the self-consistent cou- 
pling between p-h and p-p parts of the HFB Hamiltonian. 
That is, the increased probability of finding a correlated 
pair of neutrons at large distances impacts the probabil- 
ity of finding a single neutron in the halo region. It is 
to be noted that the 'halo' effect seen for q;=1/6 solely 
results from pairing, and it is not related to reduced bind- 
ing energy. In contrast, as it is shown below, increased 
pairing correlations lead to greater separation energies 
and lower chemical potentials, i.e., to increased particle 
stability. 

As emphasized in Ref. Q, pairing correlations can be 
enhanced in weakly bound nuclei due to increased surface 
effects and the closeness of the particle continuum. In 
turn, pairing can influence quite dramatically the asymp- 
totic properties of density distributions in drip-line sys- 
tems. This is nicely illustrated in Fig. ^ which com- 
pares the HFB-|-SLy4 densities calculated with a=l/2 
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for i20Sn (well bound), ^^°Sn (weakly bound), and i™Sn 
(very weakly bound) drip-line nuclei. One can see that 
adding neutrons results in a simultaneous increase of den- 
sity both in the nuclear interior and in the surface region. 
At very large distances, the asymptotic behavior of Pn{r) 
reflects the gradual rise of the Fermi level A with a neu- 
tron number. However, this effect is much stronger for 
the pairing density Indeed, as seen in Fig. ||, as com- 
pared to ^^°Sn, there is a dramatic increase in /5„(r) in 
the outer regions of weakly bound nuclei '^^"Sn and, in 
particular, ^^"Sn. These calculations indicate that for 
small values of a the box size should be chosen as very 
large if one aims for a very accurate description of HFB 
densities at large distances. 

The structure of HFB densities determines the behav- 
ior of the self-consistent p-h and p-p potentials. Figure ^ 
shows the behavior of U{r) = T{r, r) (||) (local part only 
- see discussion in Ref. [Q) and U{r) = h{r,r) (0) ob- 
tained for neutrons in the weakly bound nucleus Sn 
in the HFB+SLy4 model for several values of a. The 
behavior of the pairing potential U{r) is consistent with 
the pattern shown in Fig. |^. Indeed, for DDDI, the pair- 
ing potential is proportional to the product of the pairing 
density p and the pairing strength factor /pair (^. Con- 
sequently, U{r) is essentially peaked around the nuclear 
surface, and both its minimum and range shift towards 
larger values of r with decreasing a. For a=l/6, the 
pairing potential is still sizeable at large distances reach- 
ing 14 fm (i.e., twice the nuclear radius). The central 
p-h potential U (r) very weakly depends on the form of 
the pairing interaction. It is only for small values of a 
that, due to a direct contribution from the pairing den- 
sity to the p-h potential [see Eq. (A. 5a) in Rcf. [|l9|], a 
small barrier develops just beyond the nuclear surface. 
That is, the central neutron potential becomes slightly 
repulsive at r~9fm. However, this effect is more than 
compensated by the increased pairing potential and the 
total binding energy decreases. 

Figure H compares the HFB-|-SLy4 HFB potentials cal- 
culated with a=l/2 for I20,i50,i70gj^ With increasing 
neutron number, the radius of the p-h potential increases, 
and the potential becomes more wide in the outer region 
(i.e., it becomes more diffused). The p-p potential be- 
comes more surface-peaked and its range increases. By 
analyzing Figs. || and one can conclude that it is in 
TV-rich weakly bound nuclei that the density dependence 
of the pairing interaction is most important. 



IV. PAIRING GAPS, SEPARATION ENERGIES, 
AND HALOS 

A key feature of the nucleonic density p{r) is its second 
radial moment, i.e., the rms radius i?rms=(''^)^^^- We 
consider here and in the following the scaled rms radius. 



R 



gcom 



5 // (Pr r'^pir) 
3 V Jd^rpir) ' 



(10) 



which we call the geometric radius |24|. The factor 
(5/3)^/^ serves to bring that value closer to the box ra- 
dius. Further density characteristics are best deduced 
from the corresponding form factor 



e"''' p{r)d^r. 



(11) 



For the spherical density distribution p{r), the form fac- 
tor is also spherical and can be expressed in the standard 
way as F{q) = J jo{qr)p{r)r^dr. 

Good approximation to typical nuclear form factors is 
given by the Helm model (25|-|3^, where nucleonic den- 
sity is approximated by a convolution of a sharp-surface 
density of radius _Ro with a Gaussian smoothing profile. 



i.e. 



p(H)(r)^ / d'r'poeiRo-\r' 



exp 



f jr-r') 
{ 2^ 



(27r)3/2a3 



(12) 



The quantity Rq in Eq. (|l^ ) is called the diffraction (box 
equivalent) radius, and the folding width a represents 
the surface thickness. The diffraction radius, Rq, can be 
deduced from the first zero, qi, of the microscopic form 
factor F{q): 



Ro = 4.49341/gi, 



(13) 



while the surface thickness parameter, a, can be com- 
puted by comparing the values of microscopic and of 
Helm form factors, F{qm) and F^^^ (9m), at the first max- 
imum q,n of F{q): 



q^ RoqmF{qm) 



(14) 



The geometric radius of the Helm model can be easily 
computed as 



(15a) 



From this relation one sees that the geometric radius be- 
comes the box-equivalent radius in the limit of a small 
surface thickness. The Helm model follows the exact den- 
sity distribution over a wide range of densities, but some 
deviations may build up far outside the nucleus at very 
low densities. Thus the Helm radius is a good approxi- 
mation to the true geometric radius in well-bound nuclei. 
The situation changes if one goes towards the drip line 
where the nucleons become less bound and the outer tail 
of the density makes a non-negligible contribution to the 
geometric radius. Since the outer tail is not contained at 
all in the Helm model, the halo size can be characterized 
by the difference of these two radii |24|, i.e., 
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<5^halo = ^goom — -RHclm- (16) 

In Ref. it has been shown that (5i?haio is small in well- 
bound nuclei, but it becomes enhanced for heavy exotic 
systems with low neutron separation energies. Further- 
more, it was noticed that the halo parameter obtained in 
the HFB-|-SLy4 model with DDDI is significantly larger 
than that predicted in the HFB-|-SkP model, and in rel- 
ativistic mean-field models with the finite-range Gogny 
pairing. In the following, we shall investigate the influ- 
ence of pairing on (5_Rhaio- 

The bottom portion of Fig. || shows the average neu- 
tron pairing gaps 

(A„) = I d^rd^r'Y,hn{ra,r'a')pn{r'a\ra). 

(17) 

For a=l/6, pairing correlations are so strong that they 
give rise to a non-zero static pairing in the magic nucleus 
^^^Sn. For larger values of a, pairing gaps at 50<7V<100 
are almost independent of a, while a weak dependence is 
seen only near the drip line. 

The impact of the pairing Hamiltonian on the two- 
neutron separation energies in the Sn isotopes is illus- 
trated in Fig. ^, middle portion. The most striking result 
is that the strong pairing at low densities dramatically 
reduces the shell efi'ect at N=82. The presence of this ef- 
fect gives a strong argument against taking small values 
of a when aiming at realistic calculations. Another con- 
sequence of enhanced pairing is the shift in the position 
of the two-neutron drip line, S2n=0- While for a=l and 
1/2 the last two-neutron bound Sn isotope is calculated 
at A^=120; at a=l/6 the nuclear binding is increased by 
pairing and the drip line shifts to N—126. 

Finally, the top panel of Fig. ^ displays the neutron 
halo parameter ( |l^ ) . As already pointed out in Ref. , 
this quantity is very small in well-bound nuclei and in- 
creases with decreased separation energy. In the case of 
Sn isotopes considered, there is a gradual rise of (5i?haio 
for A'^>132. However, halo size is a very sensitive function 
of the pairing interaction. Indeed, when decreasing the 
pairing exponent from a=l to a=l/6, the halo param- 
eter increases by an order of magnitude. (Actually, for 
Q!=l/6, (5i?haio is nonzero also for well-bound systems.) 
This is a consequence of a strong feedback between p-h 
and p-p densities (cf. discussion around Fig. 0). 



The main conclusion of our work is that, due to the 
self-consistent feedback between particle and pairing den- 
sities, the size of the neutron halo is strongly influenced 
by pairing correlations; hence by the pairing parametriza- 
tion assumed. Consequently, experimental studies of 
neutron distributions in nuclei are extremely important 
for determining the density dependence of pairing inter- 
action in nuclei. Our analysis suggests that the strong 
low-density dependence of the pairing force, simulated 
by taking very small values of a in DDDI, is unphysical. 
Therefore, predictions of very large halos in drip-line nu- 
clei, obtained with a— 1/6 [g4|, should probably be re- 
visited when more reliable information on pairing forces 
becomes available. The present experimental data are 
consistent with 1/2^q;^1. In this context, it is interest- 
ing to note that excellent fits to the data were obtained in 
Refs. |l5| , |l6|] by taking a— 2/3. However, at present there 
is no theoretical argument why the density dependence 
should be even taken in a form of the power law. Clearly, 
a more precise determination of the density dependence 
(including the isovector dependence) is a challenging task 
for future theoretical work. 
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V. CONCLUSIONS 



This work contains the theoretical analysis of particle 
and pairing densities in neutron-rich nuclei and their de- 
pendence on the choice of pairing interaction. The main 
goal was to see the effect of pairing on spatial character- 
istics of nucleonic densities of nuclei far from stability, 
where the closeness of the particle continuum qualita- 
tively changes the physical situation. 
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FIG. 1. Radial strength factor /pair of the den- 
sity-dependent delta interaction, Eq. (|^, as a function of p 
for several values of a. The value of po was assumed to be 
0.16 fm~^. At each value of q, the strength Vb was adjusted to 
reproduce the neutron pairing gap in ^^"Sn. The inset shows 
/pair/ 1 Vol as a function of dimensionless normalized density 
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FIG. 2. Self-consistent spherical HFB+SLy4 local densities 
p(r) (top) and p(r) (bottom) for neutrons in ««Sn for several 
values of a. The insets show the same data in logarithmic 
scale. 
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FIG. 3. Self-consistent spherical HFB+SLy4 local densities 



p{r) (top) and p(r) (bottom) for neutrons in I20>i50,i70gjj 
a=l/2. The insets show the same data in logarithmic scale. 
The dramatic fall-off of densities at 30 fm is due to the box 
boundary conditions. 
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FIG. 4. Self-consistent spherical HFB-|-SLy4 local poten- 
tials U{r) and U{r) (shown in the inset) for neutrons in ^^"Sn 
for several values of ct. 



> 














^^^^^ ''''' 


— i2»Sn 






— i5»Sn 




iv 


--- "«Sn 









5 10 15 20 25 30 
r(fm) 



HFB+SLy4] 



5 10 15 20 25 30 

r(fm) 

FIG. 5. Self-consistent spherical HFB-|-SLy4 local poten- 
tials U{r) and U{r) (shown in the inset) for neutrons in 

120,150,170 g^^^^ a=l/2. 
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FIG. 6. Neutron halo parameters (top), two-neutron sep- 
aration energies (middle), and average neutron pairing gaps 
(0) calculated in the HFB-|-SLy4 model with different den- 
sity-dependent pairing interactions (^). 
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